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Nesushchaya  poverkhnost'  v nestatsionarnom  potoke  vblizi  _ekrana 


DATE 


LIFTING  SURFACE  IN  AN  UNSTEADY  FLOW  NEAR  A SCREEN 

[Belousov,  B.  N. , A.  N.  Lukashenko,  and  A.  N.  Panchenkov,  Nesushchaya  poverkhnost' 
v nestatsionarnom  potoke  vblizi  ekrana-.  In:  Aircraft  Construction  and  Technology 
of  the  Air  Fleet  (Samoletostroyeniye  i Tekhnika  Vozdushnogo  Flota) , No.  18, 
Khar'kov  State  University  Publishing  House,  Khar'kov,  1970,  pp.  3-11;  Russian] 

^^Several  major  works  dealing  with  the  theory  of  a lifting  surface  in  an  /'. 

unsteady  flow  with  various  approximate  methods  of  wing  design  are  now  in  exis- 
tence . It  is  common  knowledge  that  the  problem  of  a thin  lifting  surface 

of  arbitrary  aspect  ratio  is  reduced  to  two-dimensional  integral  equations  for 
which  there  are  no  closed  solutions. 

vj^0ne  of  the  important  aspects  is  the  study  of  two-dimensional  integral  equa- 
tions using  one-dimensional  approximations,  and  obtaining  final  results 
for  the  hydroaerodynamic  characteristics  of  a wing  in  an  unsteady  flow  on  the 
basis  of  these  equations. 

While  the  Prandtl  lifting  line  theory  defines  the  Prandtl  equation  as  a 
one-dimensional  approximation  for  a wing  in  a steady  fl-ow,  the  complexity  of 
the  physical  phenomenon  gives  no  unique  approximation  method  for  a wing  in  an 
unsteady  flow.  This  accounts  for  the  large  number  of  treatments  (up  to  20, 
noted  by  Bisplinghoff , Ashley,  and  HalfmanI)  which  yield  various  one-dimensional 
equations. 

Panchenkov6*7>8»9  has  developed  a method  of  integral  operators  of  the  theory 
of  a lifting  surface  in  an  unsteady  flow. 


The  general  solution  of  the  problem  is  in  the  form  of  three  components,  two 
of  which  are  obtained  from  equations  of  the  same  form  as  those  of  the  theory  of 
a wing  in  a steady  flow. 

Reference  8 also  discusses  a method  of  formulating  the  Prandtl  equation 
by  using  the  example  of  a wing  in  an  unbounded  fluid  flow.  A considerable  ad- 
vantage of  this  method  consists  in  the  fact  that  the  Prandtl  equation  was  con- 
structed for  the  parametric  constant  of  the  singular  solution.  This  permitted 
a deeper  analysis  of  the  problem  and  a critical  evaluation  of  the  existing 
theories  of  Reisner,!  Kiissner ,2  and  others. 


In  the  present  paper,  the  method  of  integral  operators  is  applied  to  the 
problem  of  a lifting  surface  near  a screen.  Two-dimensional  integral  equations 
and  a one-dimensional  integrodif ferential  equation  for  the  parametric  constant 
were  obtained.  It  should  be  noted  that  the  dovnwash  is  related  only  to  the 
singular  solution  corresponding  to  the  singularity  in  the  entering  wing  edge. 
Transposition  of  this  result  to  the  steady  motion  theory  shows  that  the  classical 
lifting  line  theory  is  strictly  applicable  only  to  a plane  lifting  surface. 

If  the  lifting  surface  is  not  planar,  terms  of  order  j^-,  not  considered  in 

the  classical  theory,  appear  in  the  expression  for  the  wing  lift  coefficient. 

From  a physical  standpoint,  they  determine  the  effect  of  the  aspect  ratio  on 
the  zero-lift  angle  of  the  wing. 


Numbers  ■'n  the  r-*gbt  margin  indicate  pagination  in  the  original  text. 
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The  present  study  is  closely  related  to  those  of  Refs.  6,  7,  8,  and  9, 
and  we  therefore  attempted  to  mention  their  general  results  very  briefly  where 
necessary,  while  discussing  the  new  results  more  extensively. 

1.  Let  a thin  lifting  surface  S move  at  velocity  Vq  near  a solid  screen.  /4 
In  addition  to  the  basic  translational  motion,  it  also  executes  harmonic  oscilla- 
tions of  frequency  oo. 

Assuming  that  the  disturbances  introduced  into  the  flow  by  the  lifting  sur- 
face are  small,  the  problem  can  be  linearized  in  a known  manner^ in  the  coor- 
dinate system  oxyz  moving  at  velocity  Vq  (see  figure). 


In  the  case  of  harmonic  oscillations,  for  the  hydrodynamic  potentials  of 
velocities  and  accelerations  0 and  their  derivatives,  we  have 


• V'igi  0 =e; (£)*“'. 

0 = o".  (g)  &mt g 6 20. 


V-  ; 
<0» 


The  index  n denotes  an  nth  order  derivative  of  the  functions  0 and  <p  with  respect 
to  the  xi  coordinate  (x3  = x;  x2  “ y;  x3  = z) . 


in 


Only  the  real  parts  of  the  corresponding  expressions  are  considered  above 
1-qs.  (1.1),  but  for  simplicity,  the  symbol  Re.  is  omitted.-  Relations  (1.1) 


' it  possible  to  formulate  the  problem  for  0(g)  and  <Hg)  and  to  exclude  time 
!>oo  consideration. 


Changing  to  dimensionless  space  ft  aL  an  incident  flow  velocity  of  1,  and 
tv  uving  the  upper  dash  from  the  corresponding  symbols  0(g)  and  <Mg) , we  have 
n boundary  value  problem  in  the  space  of  acceleration  potential  0: 


TCWWWyL'A-r  i.  zlX g.l.1 


Problem  (A)  is  given  by  the  integral  operator  Ay,  specified  in  space  L. (S  ), 
2 . 1 P 

with  the  values  C ft.  The  necessary  properties  of  operator  AY  are  given  in 

Refs.  5-8. 

In  the  problem  under  consideration,  operator  Ay  has  the  form 


where 


i j 7 (P)  i [t  + si£n  F ^dS; 

SP  .. 

r = V(x-S)»  + (g-ti}»  + (z--C)»; 
ri  = / (x  - ')'■  + {y  — v)*  + (z  + t + 4/7)»; 


sign  F. = { 


-1  - solid  screen 
+1  - free  surface. 


From  condition  (1.2),  we  obtain  a two-dimensional  integral  equation  of  the 
problem 

4*7  = f(g),  g*Sp-  ^1,5) 

The  solution  will  be  represented  in  the  form  of  three  components:^-® 

7(g)  = 7i(&)  + 12(g)  +7i(g)  g*Sp> 

where  y^  is  a solution  of  class  C^S  ),  related  to  the  presence  of  a discon- 
tinuity of  tangential  velocities  when  surface  is  crossed; 

Y2  is  a solution  of  class  C^(Sp)  describing  inertial  motions; 
y^  is  a singular  solution. 

An  analysis  of  Eq.  (1.5)  and  a method  of  setting  up  the  equations  for  y.. 
are  given  in  Ref.  6. 


We  will  write  the  equations  for  y^  and  y^: 

tf.A7»=fi'  + C.  g€S„; 

WoAlj  = —ikFi,  gsSp. 


' (1,6) 
(1.7) 


Here  F^  = - — F2(g);  ^(g)  a function  describing  the  shape  of  the  surface  S; 
N01  " ~ / dT* 

Constant  C is  obtained  from  the  solvability  condition  of  Eq.  (1.6)  in  space 

cl(V-6  • 

— 10 

Performing  the  calculations,  we.  have  for  N0.A  y 

01  z , 

1 * *H  - - _ 

v I I /(r-;)a  -|-  r-triXy-T,)^ 

-I  - I 

(r-DKv-T.inM’,'  (71  ->’  <i>  I6A5 Wl  (’|)  I*5** Will  ,c 

—sing  h } 

'x  I 

. (1,8) 


; 


where 

TW-ift- 

The  singular  solution  is  determined  from  the  equation 
NAj  = Fx  (gY,  h € Sp\ 

X 

N = — 4kx  ( e~lk'd-,  0$) 

•* 

where  k is  the  chordwise  Strouhal  number. 

Equation  (1.9)  is  obtained  by  mapping  Eq.  (1.5)  into  the  space  of  the  /6 

velocity  potential. 

2.  Let  us  now  consider  the  Prandtl  problem  for  a wing  in  an  unsteady  flow. 

For  a wing  of  high  aspect  ratio,  the  operator  N01A2y  is  approximated  by  the 
operator 

. ‘ 4-1  ' “it  ' 

_ ■ i/Jai ••  r . . ..  ■»  ivaiw* 


V‘T  ^ | T ft  -»  [f=6  + .*„-£+ W>  J |f| 


The  form  of  singular  solution  y^  is  known  for  plane  flow,  and  "for  a wing 
of  high  aspect  ratio  we  will  take  y3  in  the  form  of  a solution  to  the  two- 
dimensional  problem,  but  with  a constant  dependent  on-  the  y coordinate: 

. 7* ~a(i/) y j±$. ,/ . ..  •••  . . 


we  have 


l°l2  6 C1  (S)  h 7,  = —ik  j 7,-cfc;  g £ Sp 


na, rJa  = )v0,4ti;  S € Sp. 


Considering  Eqs.  (1.6),  (1.7),  (1.9)  and  the  property  (2.2),  we  obtain 
NAji  = — C — W J,T«- 

Equation  (1.6)  can  be  written  in  the  form 

+ O,  (2*3? 


".--J - • ii 

‘■a 

For  X -*-®C1  C,  and  for  higher  aspect  ratios,  we  will  take  the  approximatit 

+«  ‘ VV1S* 

C = —1 1*  tf»— (2.4)' 
n J pi  — 


-(V. dx 


Then,  introducing  the  operator  p : 


>",-Vo,/./,»  - { Ti  (*)<?' 


and  taking  approximation  (2.1)  from  (2.3),  we  obtain  the  Prandtl  equation  for 
a regular  solution  having  the  same  form  as  the  Prandtl  equation  of  the  stationary 
problem  :*•  7 


where 


’lit 
' “j 

• +1 

A = j T.  ($)<*;.  '■  S 


The  function  Y22  obtained  from  the  equation 

_ 1 

= ikC(y).  ‘ (2.6) 

We  will  introduce  a series  of  assumptions  typical  of  the  Prandtl  problem:^’8  /7 

~ Nl0Aj  + NnAj; 

7 (P)  = 7i  (?)T2 (y);_  ' : " 

aT.  = 27*  (y)  >•  (y)  A^jo  A,ti  («); 

• t 

iVu  -=  — c1** ^ e~tMxtk;  ..  .... 


Ar,0<4,7  (i  7 (')  L - « + Sign  ^ 16S*(J  ^ + 


— I 

+i  x 


■ ,+ike‘^i  7 (')  J e~‘k'  [(7~r|j  + signF— 


(*~i)'+lW(y)\ 


«(y)2>  (y)  A(10/l,Tl  (')  = -C  (y)  - 21  (y)  Tm  (y)  Nl0A2  Tin  - 

-A/i^,[u(i/)]/i±|+7J,];  £€SP.  (2.7) 

F°r  ikxln8  ln  ^ Unbounded  fluid.  b°th  parts  of  Eq.  (2.7)  contain  the 

factors  e . Multiplying  Eq.  (2.7)  by  e~lkx,  we  eliminate  the  variable  x.5,8 
This  property  also  applies  in  the  general  case,  but  the  equation  cannot  be 
ilex 

solved  for  e in  explicit  form,  since  one  can  calculate  the  integrals  in 

N10Az7  in  closed  form.  We  therefore  introduce  the  averaging  of  Eq.  (2.7)  in 
accordance  with  the  rule 


— I 


^or  and  7«*'(y) /V,„]4^ni  VG  will  use 


N\sA,-\x  « 


«k,  (*.  2* 


T.„  M Af,  JL.  - ^ |t;  ■ ; [| ».  we  m + C,  MU,  w]  ■ (2.9) 

For  h -*■  °°,  the  representations  (-2.8)  and  (2.9)  change  into  the  expressions 
obtained  in  Ref.  9.  The  functions  N^k)  were  also  obtained  in  Ref.  9. 

Then  t,  ... 


Nx(» 

. 7jii  (y)  i 

F_£*F  A>  (k)rtrA 

•;,(*;  h)  2ir 

r 2«J<,  (A:  h)  '("2j  (y)  | 

= - 

-C(j/)-JVu"/4,[a  ([/)  | 

/TT5  , 

/ 7«| 

+i  

Upon  introduction  of  the  new  variable  (!/)  = J !< 2 (y)  v\  + 7221  (y;  5)1  dh  , / 

— i ■ ■* . 

E«j.  (2.8)  becomes 

. MO 

' ^iW  ; 

+1  . - * ^ 

<W  (f) « ^ 1 ^ 

. . , -i  . ■■ 

As  follows  from  the  solution  of  the  two-dimensional  problem  of  unsteady 
motion  of  a wing  near ‘a  screen, /»B  the  function  y(S)  may  be  written  as 

7 (S)  = ^a0(y)  ]/i±|  + *7?  (5)  + (S). 

7ne  functions  Y^(S),  corresponding  to  the  two-dimensional  problem  for  an  un- 
bounded fluid,  were  calculated  in  Ref.  9.  Using  these  results,  we  have 

[-  ( k ) C (y)  -f  Cj  (y)  Ni  (ft)]  — l)  ik C3  (y)^,  (2,12)  j 

LaX  being  an  integral  operator  whose  kernel  was  obtained  in  the  monograph  of  Ref.  4. 

Va  must  now  establish  the  relations  between  the  functions  excluding  and 
v2  ^loin  Kq.  (2.9).  By  definition,  when  X -*•  « 


„ /*,«#»>,  +,OrC,(y)l 

“•m— |-c«r+~nsrj- 


The  constants  C^(y)  are 


c.v>- H 


t-'  

O <»--■;$  V TTifMd‘- 


From  Eq.  (2.9)  we  find 


aI  = “j.  = + i*‘* C3  (y)  — 

X--  Yf 

- [£  Rt  (k)  C(y)  + C3  (y)  N,  <*)]  (|  - 1) . 

Comparing  the  two  expressions  for  a^(y),  we  obtain  the  desired  relation 

fy“4W  ^i  = V 

Combining  the  results  obtained  from  Eq.  (2.9),  we  write  the  integrodif fer-  /9 

ential  equation  of  a high-aspect-ratio  wing  near  a screen  in  an  unsteady  flow 

- +1 

M - x|i [». - -snifcwi  1 “!■  fo  K Ifa (#  - >1)1  *1.  (2,13) 


where 


The  kernel  K (9)  is4 


» — 'i 

(</  - V + 16** 


- n_i[p  V(y—W  +.i6ft»- 


(y  — ri)!  + 16** 


No[pV  (y-rly-+  16/i2). 


Wljpn  K 0,  N^(k)  -*■  -1,  and  from  Eq.  (2.11)  we  obtain  the  Prandtl  equation  for 
a wing  near  a screen  in  a steady  flow.4  Then 


K0/-^)  = (7^)  + signf^ 


(y—,) 

■ T,)’  + 16**  ■ 


Thus,  the  Prandtl  problem  for  a wing  in  an  unsteady  flow  has  been  reduced 
to  two  one-dimensional  integral  equations  (2.5)  and  (2.13),  with  Eq.  (2.5)  cor- 
responding to  the  equation  of  the  stationary  theory. 

This  result,  which  is  a direct  consequence  of  the  separation  of  the  general 
solution  into  three  components,  is  of  major  importance,  since  by  introducing 
"reduced  circulation"  type  integral  quantities,  the  authors  of  the  known 
theories  find  a single  one-dimensional  equation.  v 

3.  In  determining  the  aerodynamic  coefficients  C^  and  C^,  a close  first 
approximation  is  given  by  the  solutions  given  in  Refs.  3 and  4. 


r (y)  y\ 


We  will  use  this  representation  for  an  approximate  solution  of  (2.13) 
a . (y)  = <I>  1/  1 — y1  we  obtain  for  constant  $ 


Using 


2X(0)lN,(i)l 


where  £k^  is  given  by  the  formula  from  Ref.  4: 

+i  +i 

«>-  7T=f  -5=9)  «V<W> 

k^  being  the  spanwise* Strouhal  number. 

The  wing  lift  coefficient  is  calculated  from  the  formula 


-t-i  -t-i 

Cv  *=4  J J 7 ('•  Ti)d\dn. 


We  will  write  the  wing  lift  coefficient  in  the  form  of  two  components: 

Cy  — Cyv  -f-  Cvi  > 

where  is  the  lift  coefficient  determined  by  the  vortex  motion; 

is  the  lift  coefficient  determined  by  the  apparent  wing  mass. 

We  find  the  components  C and  C 

yv  +i  yi 


■ r — ' ' C(«  f ■}/  1 — x C0  , \ J . 

- ~rs^r 

^ •:  + *x|iv,(*)  | c*>. 

;+ * [ 7T^xdx\'-  tt— f it ; 


c,«  - 2 ik \rTZTfFl(x) dx.  (sS 

Here  ' ''  ; ' ■ • : • : ^ 

. . ; Z (*)  - F,  (x)  ~ [A'0  - AUATi.  IS 

ry  • 

C(k)  being  the  Theodorsen  coefficient. 

in  the  calculations,  the  theoretical  value  a<u  = 2tt  (allowing  for  the  in- 
fluence of  the  viscosity  of  the  fluid)  may  be  taken  as  aro  ~ 5.45. 

For  sCeadv  flow  C , » 0,  and 

yi  . 


8 


Cj,-V  


ij  ] V^/-^+ 


+ 2>  f — xdx  1 r -1 . 


The  first  term  in  Eq.  (3.6)  gives  a known  value  of  the  wing  coefficient 

in  e bounded  flow,  and  the  second  term,  of  order  is  the  refinement  introduced 

into  the  Prandtl  theory  by  the  present  method.  For  a plane  lifting  surface,  the 
second  term  in  Eq.  (3.6)  is  zero,  and  the  results  will  fully  correspond  to  the 
classical  theory  of  a lifting  line.3>4 

The  functions  \p^,  ip  were  obtained  in  Refs.  7,  8,  and  11  for  translational 
and  rotational  wing  oscillations. 
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